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APPLICATION OF CONE-GRID SCHEME FOR SOLVING BURGERS 
EQUATION 


MAHMOUD A.E. ABDELRAHMAN 


ABSTRACT. The cone-grid scheme for the Burgers equation is introduced, which gives exact 
unique solutions. The Riemann solution for the Burgers equation is the building blocks for 
this scheme. To validate this method and assess its accuracy, front tracking scheme is 
also applied to the same model. These two schemes are compared by interesting numerical 
examples, where the explicit solutions are known. 


1. INTRODUCTION 


In this paper we consider the Burgers equation, which given in the following form : 


(1.1) Uz + E = 0, 


which we can write in the form 
(1.2) Ut + Uug = 0: 


The numerical solution of Burgers equation is so importance due to the equations appli- 
cation in the approximate theory of flow through a shock wave travelling in a viscous fluid 
and in the Burgers model of turbulence. 

The front tracking and the cone-grid schemes could serve as a common numerical technique 
to solve several scientific and engineering problems due to avoiding the specific eigenstruc- 
ture of the problem [1]. In fact the Riemann solution for the Burgers equation is the basic 
ingredients for these schemes. These schemes have also high order accuracy in smooth re- 
gions of the simulated flow and resolve sharp discontinuous profiles in the shock regions. 


In [1] we presented a new schemes for the first time, namely the cone-grid scheme. Here 
we will show that the Burgers equation equations give an surprisingly explicit algebraic solu- 
tion. The cone-grid scheme also gives a powerful shock resolution. Hence the cone-grid and 
front tracking schemes enable us to give an instructive comparison with each other as well as 
with well known exact solutions, for example the exact Riemann solutions. Two numerical 


examples studies are carried out. 
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In this paper, the cone-grid scheme is implemented for solving Burgers equation in one 
space dimension. For validation and comparison, the numerical results of front tracking 
schemes is also presented [3]. 

The paper is organized as follows: Section 2 gives a brief description of the Burgers 
equation. The one-dimensional cone-grid scheme is briefly presented in Section 3. In Section 
4, one-dimensional numerical test problems are given. Finally, in Section 5 the conclusions 
are summarized. 


2. THE RIEMANN PROBLEM 


A Riemann problem for the Burgers equation, consists of a conservation law, 


(2.1) u + E = 0, 


together with a piecewise constant data having a single discontinuity, 


=, <0 
(2.2) utr) ee 
u, Wed: 


The form of the solution depends on the relation between u_ and u4. 


e Casel. If u_ > u4, in this case there is a unique weak solution, 


_ E< 
(2.3) DCH EE ES, 
Ur T >S 
where 
(2.4) ee 


2 
is the shock speed, the speed at which the discontinuity travels. 

e Case 2. If u_ < u‚, in this case there are infinitely many weak solutions. One of 
these is again (2.3) and (2.4) in which the discontinuity propagates with speed s. 
Another weak solution is the rarefaction wave 


u, $ <u 
(2.5) u(x,t) = Z u <Ẹ<uz , 
Ur, T > uz 


For more details about the Burgers equation we refer to the textbooks [2], [4] and [5]. 


3. CONE-GRID SCHEME 


In this section, the cone-grid scheme is derived [1]. The one-dimensional Burgers equation 


for numerical purposes may be written down in the dimensionless “vector form ”, 


(3.1) ui + Fil, = 0. 
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FIGURE 


FIGURE 1. Balance re- 2. Computational 
gions [1] . domain [1]. 

In this section we review the cone-grid scheme, which is discretized in time and space. This 
scheme is based on the integral conservation laws in terms of curve integrals acclimatized to 
the choice of the numerical time and space grids. Moreover we will prove that this scheme 
admits a unique weak solution for the Burgers equation. 


This scheme is very interesting and can be obtained easily from the conservation laws on 
a cone-grid. We describe a fixed time step At > 0 and calculate the spatial mesh size Ax in 
terms of the natural Courant-Friedrichs-Lewy (CFL) condition, which is satisfied. 

We consider the Riemann solution of the Burgers equation inside the cone depicted in 
the Figure 1 in order to derive the cone-grid scheme. For the construction of the cone-grid 
scheme, we use the conservation laws with respect to the computational domain depicted in 
Figure 1, namely 


(3.2) puc, x)dx — F(u(t,x))dt = 0. 


on 


3.1. The computational domain. Before applying the scheme, it is necessary to discretize 
the computational domain. The computational domain is the trapezium in Figure 2. Recall 


AS max fe) Ie 


a<u<b 














Given are a and 8 € R, with a < 8, T > 0 and M EN, which give the number of grid cells 
for the half of the space interval [a, 6], i.e. 


(a) The spatial mesh size is 


(3.3) Art := 
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FIGURE 3. Balance cell|1]. 


(b) The number 


zo NE 


is the smallest integer such that ÀA. T < N - Az, where the ceiling function [y] is the 





smallest integer larger than or equal y. 
(c) Put At := 5. Then T = 2NAt. Hence, using 1 and 2 we have 


(3.5) ALL = 





where the CFL-condition is satisfied. N gives the grid cells number’s for the half of 
the time interval [0,7] and At the time step size i.e. the interval [0,7] divided into 
2N sub intervals of length At. 

(d) We have the time discretization t; = j- At for j € {0,...,2N}. 

(e) We define the quantities 


(3.6) vW=a-T-A, 8=8+T-A, 


which give the bounds of the trapezium domain on the z-axis to avoid the artificial 
numerical boundary effects, see Figure 2. The spatial intervals Jo, ol, 10. 8'] split 
into N sub intervals of length Az. 

(f) At time t;_1 we define the grid points 


£j =@ +(j-1)A-At+(k-1)At, j=1,...,2N4+1, kel, MAN) -—j+2. 


We note that the quantities x; determine all points of each balance cell in Figure 
1. Especially we can easily check that zj} = a’ and aa = Ó at initial time 
to = 0 as well as zon+11 = Q and Toni am41 = H at final time ten = T, which give 
the four corner points in the trapezoidal computational domain 2. 
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3.2. The construction of the solution. In this section we consider a restriction to spatial 
grid points that gives the left corner points of each balance cell in Figure 1 as well as in Figure 
3. In these balance cells we will solve the local Riemann problems. For the left corner-point 
vı given by vı = (tj—-1, rat) in Figure 3 we solve the local Riemann problem 


ut + Uur = 0, 


eS E e 
u(tj—ı, £x) = 

u} if HDE 
for t > t;-1 and assume that the numerical solution inside the parallelogram is given by 
u(t, x). 
Especially the constant value of u(t, x) along the cord connecting vı with v3 will be denoted by 
u. Consider the triangular balance regions OF, and Q; , spanned by the points vi, 13, va 
and vi, vo, V3, respectively, see Figure 3. On each cord 12v3 and sat we replace the 
numerical solution by the unknown constant values u_ and u, which require that 





(3.7) J PoE 
d 
For the balance region ORT, using 3.5 we obtain: 
/ 1 / 
(3.8) Uu [2 F(i) — F(u’,) — F(u,)| =0, 





d . D 1 
i.e. an implicit equation for u. 


For the balance region 09; ;, we obtain in the same way the following implicit equation for 


U: 


(3.9) u — u + R [2 F(a) — F(u 


The new states u- turn out to be uniquely calculated by using 3.8, 3.9 and the CFL-condition 
Ae all 





Note that for fixed values of a, 8 and T the quantity A in (3.5) depends on M, such that 
we can also rewrite it in the form A = Am. Then we conclude from (3.3), (3.4) that 


fm Awad 
Complete derivation of the scheme introduced in [1]. 


In the next Lemma we show that the solution given by the cone-grid scheme is unique. In 
fact we don’t prove this interesting result for the ultra-relativistic Euler equations [1]. 
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Lemma 3.1. Let be u_ < ur. Then equation (3.8) has a unique solution Uy € [ü, u4], and 
(3.9) has a unique solution u_ € [u_, ù]. 


~ 1 


Proof. For o = u, it follows u- = u, = D = u} u_, hence assume u_ < u,. By the 








formal replacements F > —F, u_ > u} and u_ > u, equation (3.9) reduces to (3.8). 
It is thus sufficient to prove the existence and uniqueness for Uy :Assume first that u, < u, 
are two different solutions. Then 











H / At U / 
Us Uy = arr (u) — F(u,)) >0, 
and 
Wi 1 At 1 / 
up-u, = ~ 25 E (+) — F(u) 
SC At F(u,) = F(u.) Hi b 
E SE Uy u, (u+ = u4) 
At Hi 1 

< Inr Aged) 
< (wy, — U4) 


from the CFL-condition, which is a contradiction. It remains to show the existence for Uy : 
It follows from the convex-hull construction of ú that 




















(3.10) F(ù) = Pe F(u). 
Define B: [ü,u‚] — R by 
Als 
Blu) =u- uy — 2757 E ®© — Flu) — Flur) 
Then from the CFL-condition: 
R P At 7 
Ba) = Au +22 (Fu) - FW) 
S At id 
< Ü—-ur +2 A Ag ut — f) 
A 
mis ni en 
Ag 
>0 
(3.11) 
and from (3.10) : 
At ~ 
(3.12) Blu.) = 42 (Flus) — FW) > 0. 


Then from (3.11), (3.12) and the mean value theorem we obtain Dia, = 0 for an appropriate 











u, € [ū, u4]. 
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solution u (cone-grid scheme), N = 400 time steps 














01 02 03 04 05 06 OF 08 09 
time t 


solution u (front tracking scheme), Nt = 400 time steps 














01 02 03 04 05 Dë OF 08 09 
time t 


u(1,x), Nt = 400 time steps 
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Exact Riemann solution 
——-—Cone Grid Scheme 
—-+--Front Tracking Scheme 




















FIGURE 4. A single shock solution. 


8 MAHMOUD A.E. ABDELRAHMAN 





t=1 


FIGURE 5. Two different initial data give the same solution at t= 1. 


4. NUMERICAL EXAMPLES 


In this section we introduce numerical examples for the solution of Burgers equation. We 
use the exact Riemann solution, cone-grid and front tracking schemes for the comparison. 
We will show that how the cone-grid is powerful scheme to solve Burgers equation. 


Example 4.1. Single shock solution of the Burgers equation. 


In this example we test our cone-grid and front tracking schemes for a single shock problem. 
We supplied initial data to the program for which we know that a single shock solution results 
from the Rankine-Hugoniot jump conditions. We select the initial data and the space-time 
range such that the shock exactly reaches the left lower corner at the time axis. Figure (4)12 
represent the plots of the solution in the time range 0 < t < 1 and in the space range 
(zz The figures shows that cone-grid and front tracking schemes captures this shock 
in exactly the same way as predicted by the RankineHugoniot jump conditions. The Figure 
43 presents the solution at the fixed time t = 1 for the same initial data. The Riemannian 
initial data is chosen as 


u = 0.25 x (1 — sign(z)) 


where 400 mesh points are considered here. In this example we found that our cone-grid 
and front tracking schemes gives a sharp shock resolution. This is a good test for the front 
tracking scheme, and its success indicates that the entropy inequality is satisfied. 


Example 4.2. An example of the irreversibility for the Burgers equation. 


Irreversibility means the following: If we know the solution at one time, then we can 
reconstruct the solution in the past. In Smoller [5, Chapter 15, 8 E] one will find a simple 
example for the irreversibility for the Burgers equation. In this example we test our cone-grid 
scheme. 
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solution u of u,+uu,=0, N = 500 states 


position x 





o 0.5 1 15 2 
time t 


FIGURE o Burgers solution for 0 <t < 2. 


We define u,(t,x) to be a solution of the equation (2.1), which present the compression 


wave, for 0 <t < 1 we set 


1, mtd 
(4.1) mime, = t= leet, 
—1, x>1-t. 


We define uo(t,x) to be a solution of the equation (2.1), which present the single shock, 


we set 


1 x<0 
4.2 ie) = g = 
SS walt, e) —1, x>0. 


The ui. are depicted in Figure (5). It is easy to check that they are all solutions, moreover, 
each of u12 satisfies the entropy condition. 
Thus they are all correct solutions. The point that we wish to emphasize here is that all 


of these solutions coincide att = 1: 


1 x<0 
4.3 Lais : ani 
SS Me i —1, z >Q. 


Att = 1 we know that a shock has formed. We don’t know when it was formed, nor can 


we even say how it was formed. 
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We emphasize again that all of these solutions are the right ones, in that they belong to the 
class of solutions, which are uniquely determined by their initial values. But this uniqueness 
is only in the direction of increasing t. Two solutions in this class which agree at t = 1 must 
be equal for allt > 1, but need not be equal fort < 1. It is in this strong sense that 
solutions of conservation laws do not have backward uniqueness. 


The numerical simulation in Figure 6 serve for illustration the non-backward uniqueness 
using cone-grid scheme. 


5. CONCLUSIONS 


In this paper we described the cone-grid scheme for the Burgers equation numerically. 

Moreover, we have shown that this scheme admits a unique solution. The cone-grid scheme 
also gives a powerful shock resolution. 
We implemented the cone-grid and front tracking schemes numerically. Hence, we introduce 
numerical test examples for the solution of the Burgers equation. For the interesting com- 
parison we use the exact Riemann solutions for the cone-grid scheme as well as for the front 
tracking scheme. Finally, the numerical examples are given. 
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